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The spin accumulation due to the spin current induced by the perpendicular temperature gradi- 
ent (the spin Nernst effect) is studied in a two-dimensional electron system (2DES) with spin-orbit 
interaction by employing the Boltzmann equation. The considered 2DES is confined within a sym- 
metric quantum well with delta doping at the center of the well. A symmetry consideration leads to 
the spin-orbit interaction which is diagonal in the spin component perpendicular to the 2DES. As 
origins of the spin current, the skew scattering and the side jump are considered at each impurity 
on the center plane of the well. It is shown that, for repulsive impurity potentials, the spin-Nernst 
coefficient changes its sign at the impurity density where contributions from the skew scattering and 
the side jump cancel each other out. This is in contrast to the spin Hall effect in which the sign 
change of the coefficient occurs for attractive impurity potentials. 

PACS numbers: 72.25.Dc, 73.63.Hs 



I. INTRODUCTION 

The spin Hall effeclP^ is the generation of the spin ac- 
cumulation, or the difference in density between spin-up 
and spin-down electrons, due to the spin current driven 
by the perpendicular electric field. This transverse ef- 
fect is produced by spin-orbit interaction in the absence 
of magnetic field. It has attracted much attention in the 
field of spintronica^ as a promising way to create the spin 
accumulation in nonmagnetic materials. The first report 
on the observation of the spin Hall effect has been made 
by Kato et alP for three-dimensional electron systems 
(3DES) in semiconductors, n-doped GaAs and n-doped 
InGaAs, and is followed by many experimental works in- 
cluding the observation in two-dimensional hole systems 
(2DHS)P and that in two-dimensional electron systems 
(2DES)P Theoretical proposals have been made before 
such observations and are classified into the intrinsic ori- 
gin and the extrinsic one. The intrinsic spin Hall effeclp] 
is due to the spin-orbit interaction induced by the crys- 
tal potential as well as the confining potential of a quan- 
tum well. The extrinsic spin Hall effectP^HH] originates 
from electron scatterings from nonmagnetic impurities 
in the presence of the spin-orbit interaction. The spin 
Hall effect observed in the 3DElP and that in the 2DES 7 
have been explained by calculations based on the extrin- 
sic mechanismP^H^ In this paper we investigate the ex- 
trinsic spin Nernst effect in 2DES. 

The observation of the spin Hall effect in 2DES has 
been made by Sih et alP in a (110) AlGaAs quantum 
well. They have already suggested in their paper that 
the observed spin Hall effect is extrinsic, since (1) the 
quantum well is doped at the area density of 10 12 cm~ 2 , 
(2) the measured value of the Rashba coefficient is small, 
and (3) the Drcsselhaus field should be absent because of 
the current orientation along the [001] axis in the (110) 
quantum well. Since the measurement in 2DES, 7 as the 
3DES experiment^ is made at the temperature of 30K, 
the phase coherence in the electron transport may not 



be important. Therefore a theoretical study for this ex- 
periment has been performed based on the Boltzmann 
equation by Hankiewicz and Vignalep^ as well as the 
semiclassical theory by Engel et alP^l for the 3DES ex- 
periment. 

In the atomic-layer epitaxial growth of a semiconduc- 
tor heterostructure, both positively and negatively ion- 
ized impurities can be introduced at a precise distance 
from the heterointerface by employing the method of 
delta-dopingP^In f ac t ; both Si (donor) and Be (acceptor) 
have been doped successfully at a precise distance from 
the interface of a GaAs/ AlGaAs heterostructure, and a 
strong dependence on the dopant type has been found in 
magnetotransport properties of a 2DES located near the 
dopant. SI Such an accurate control of the doping profile 
gives the 2DES an advantage in that this can provide a 
method to enhance strongly the spin accumulation due 
to the extrinsic spin Hall effect. 

A remarkable dependence of the extrinsic spin Hall cur- 
rent on the impurity-limited mobility has bee n foun d in 
a model of 2DES by Hankiewicz and othersPEH The 
2DES in their model has a negligible width and therefore 
the dependence on the above-mentioned doping profile 
is beyond the scope of their works. There are two con- 
tributions to the extrinsic spin Hall current. One is the 
contribution from the skew scatterin g 20 " 22 ! anc j t ne other 
is that from the side jumpP^HUl Both have long been 
studied in the theory of the anomalous Hall effect in fer- 
romagnetic metals (see Refs.26 28 for early theories on 
the anomalous Hall effect and Rcf . 29 for a recent review) . 
The skew-scattering contribution has a different sign de- 
pending on whether the impurity potential is attractive 
or repulsive, while the side-jump contribution is inde- 
pendent of both the impurity potential and the impurity 
density. For attractive impurity potentials, the contribu- 
tions from the skew scattering and the side jump are op- 
posite in sign. Therefore the direction of the spin current 
is switched as the weight of the skew-scattering contribu- 
tion is changed, for example by varying the mobility! 15 ^ 19 ! 
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This theoretical finding suggests that the spin accumula- 
tion due to the extrinsic spin Hall effect can be controlled 
in a wide range, for example, by changing the impurity 
density. We expect that the controllability should be en- 
hanced by introducing various doping profiles with the 
delta-doping technique. 

The temperature gradient is another driving force for 
the spin current in the perpendicular direction. This phe- 
nomenon, called the spin Nernst effect, is one of the most 
important subjects in 'spin caloritronics', a research field 
exploring the int erplay between the heat and the spin 
degree of freedom,™!] The spin Nernst effect is the non- 
magnetic analogue of the anomalous Nernst effect. While 
the anomalous Nernst effect has been studied in 3D fer- 
romagnetic metals for nearly a century (see RefsOU and 
l33l for early experiments, Refs|24] and [34] for early the- 
ories), studies on the spin Nernst effect have started 
quite recently. An experimental study to observe the 
spin Nernst effect is in progress in 3D metalsPSl Several 
theoretical studies on the spin Nernst effect have been 
made in 2DES PaESI However, these theories are only for 
the intrinsic origin due to the Rashba term. The spin 
Nernst effect with the extrinsic origin is worth studying 
theoretically, in particular, the dependence on the type 
and the density of impurities. Even the sign of each con- 
tribution in the extrinsic mechanism is not known in the 
spin Nernst effect. 

In this paper we study theoretically the spin Nernst 
effect in 2DES based on the extrinsic mechanism by em- 
ploying the Boltzmann equation. In particular, we pro- 
pose an efficient method to control the spin Nernst effect 
by changing the impurity type and density. 

In Sec. |n]we describe our formulation. We start from 
the Hamiltonian for an electron in a quantum well formed 
in a semiconductor heterostructure with interfaces par- 
allel to the xy plane. Then we reduce it to the effec- 
tive Hamiltonian for the two-dimensional electron mo- 



tion in the xy plane (Sec. II A). Here we show that the 



2D Hamiltonian becomes diagonal in the z component of 
spin when each impurity is located on the center plane of 
a symmetric quantum well. For such 2D Hamiltonian we 
write the Boltzmann equation and derive the distribu- 
tion function (Sec. II B). Using the distribution function 
we obtain the current densities and the transport coeffi- 
cients (Sec. |II C ). We show here that the side jump also 



gives rise to the current density component induced by 
the temperature gradient. 

Then we apply the formulation to the spin Nernst effect 
in Sec. |III| We consider a rectangular 2DES, apply the 
temperature gradient along the x direction, and calculate 
the gradient along y of the chemical-potential difference 
between spin-up and spin-down electrons. We pay a spe- 
cial attention to the signs of contributions from the skew 
scattering and the side jump. We present the result as a 
function of the impurity density for both attractive and 
repulsive potentials and compare it with that of the spin 



II. FORMULATION 



2D Hamiltonian 



We consider conduction-band electron states which are 
bound to a quantum well with translational symmetry in 
the xy plane. We assume that the wave function de- 
scribing the motion along the z direction is frozen to the 
ground state, and derive the effective Hamiltonian for the 
2D motion in the xy plane in the following. 

We start from the Hamiltonian describing the 3D mo- 
tion: 



-H3D = 



Pi 



-PjL 
2m 



2 >Pl 



+ V 3D (x,y,z) - acr • (VV3D X p) , 

(1) 

where m is the effective mass, a is the effective cou- 
pling constant of the spin-orbit interaction for an elec- 
tron in the conduction band of the semiconductor, and 
c = {<?x, Cy, &z) is the Pauli spin matrix. The potential 
energy is 

V 3 D(x,y,z) = V wcn (z) + V iiap (x,y,z) + eE r, (2) 

where V wc \\{z) is the well potential, Vi mp (x, y, z) is the 
potential due to randomly-distributed impurities, E = 
(E x , E y , 0) is the in-plane electric field, and e > is the 
absolute value of the electronic charge. 

We define the Hamiltonian for two-dimensional motion 



H2D = (H; 



3D) 



(3) 



where the brackets represent the average with respect to 
the motion along z as 



(#3d) = J dz Lp (z)H 3D Lp {z). 



(4) 



Here <fo(z) is the wave function of the ground state at 
energy Eq which satisfies the Schrodinger equation: 



2m 



U W eii(z) 



ipo(z) = s ip (z). 



(5) 



We begin with evaluating terms in H2D which originate 
from the spin-orbit interaction. We here assume that 
V we n(z) is symmetric with respect to the center of the 
well, z = 0. Then V we \\(z) gives no spin-orbit term in 
i?2D- The in-plane electric field gives a spin-orbit term 
with a z only (no terms with <r x and a y ), since (p z ) = 
and E z = 0. 

Spin-orbit terms in i?2D, which is due to the impurity 
potential, are separated into the following three compo- 
nents: 



Hall effect. Conclusions are given in Sec. IV 



Trso.imp 

rrSOjmp 

U 2B,y - 

Trso.imp 

U 2D,z — 



-aa x [((VyVi mp )p z ) - (V z Vi mp ) p y ] , 

<y [(V z Vi mp ) Px ~ ((V x V imp )p z )} , (6) 

'Z [{V X V imp )Py ~ {VyV imp )p x } , 
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with the effective impurity potential in the 2DES, 



-amp 



(x,y) = (V imp (x,y,z)) 



(7) 



Since a term in -ff^D™ 13 can be rewritten as 
((V y V imp )p z ) = ih{V v (V 2 Vimp))/2 and the same is 
true for -ff2D™ P > ^ ne magnitude of #2D™ P an< ^ that 
of H 2D,y lP are determined by (V z V imp ), V x (V z Vi mp ) 
and V y (V z V[ mp ). On the other hand the magnitude of 
^2dT P * s determined by V^fimp and V y Vi mp . 

Equation ^ demonstrates that the 2D Hamiltonian 
for 2DES formed in a quantum well, in general, con- 
tains in-plane components of spin, a x and a y , due to 
the combined action of the impurity potential and the 
spin-orbit interaction. The re sulting spin relaxation due 
to the Elliott- Yafet mec hanis m 1 4 * 39 ! 40 ! has already been re- 
ported in the literature! 41 * 42 ! However, the z component 
of spin, er 2 , is conserved when the condition 



(V z V imp ) = 0, 



(8) 



is satisfied. This condition is satisfied when impurities 
are located on the center plane (z = 0) of the symmetric 
quantum well. Such a precise placement of impurit ies is 
in fact possible by using the method of delta-dopingffZEH 
We therefore assume the condition Eq.Q. Our Hamil- 
tonian for the two-dimensional motion of the 2DES is 
simplified to become 



p 2 + p 2 

#2D = ~ -+V2B(x i y)-aa z [(\' x V2D)Py-(^yV2B)Px 

2m 



with 



V2n(x,y) = v lmp {x,y) + e(E x x + E y y). 



(9) 



(10) 



This 2D Hamiltonian coincides with that employed to 
study the extrinsic spin Hall effect of 2DES in the previ- 
ous theory.QS 



B. Boltzmann equation and the distribution 
function 



Hankiewicz and Vignale in their study on the extrinsic 
spin Hall effect of 2DES 15 have obtained the distribution 
function by solving the Boltzmann equation up to the 
first order of the electric field E and of the spin-orbit 
coupling constant a. Here we extend their formulation 
to include gradients of the chemical potential and the 
electron temperature as driving forces, and obtain the 
distribution function up to the first order of all the driv- 
ing forces, which is denoted simply by 0{E) below, and 
up to 0(a). We show that the side jump, as well as the 
skew scattering, gives a temperature-gradient term in the 
distribution function. 

Since our 2D Hamiltonian conserves the z component 
of spin, the distribution function for each of its eigenval- 



ues a = ±1 is determined independently by the Boltz- 
mann equation. The Boltzmann equation for the distri- 
bution function of electrons with spin a, f a (r,k) in a 
steady state is 



dfa 
dr 



(-e)E 0% 
h ' dk 



dfa 

dt 



(11) 



The distribution function is decomposed into that in the 
local equilibrium, J*- ', which depends on k through the 
energy e k = h 2 k 2 /2m, and the deviation in the first order 

of the driving forces, which depends on the direction 
of fe relative to E: 

U(r,k) = /<°> (e fc) Mr),T e (r)) + /W(r,fe) 1 (12) 

where f °\e,n,T) = {exp[(e- fi)/k B T] + l}" 1 , (x a is the 
spin-dependent chemical potential, and T c is the electron 
temperature. Note that the first term of fair, k) includes 
spatial dependences of and T c , although the function 
/(°) itself is of the zeroth order of the driving forces. 
The r dependence of fa (r,k) also originates from the 
driving forces, and therefore it gives only terms of 0(E 2 ). 
Since 



hk 

m 



v^-+0(E), 



and 



dfa 

dr 



VT C + 0{E 2 



(13) 



(14) 



then the left hand side of the Boltzmann equation Eq.( 11 ) 
is written in the first order of the driving forces as 



dfa 

dr 



(~e)E df a 0/<°> 

-IT' 9fe = t " F " (£fe) &P 



with a generalized force 

F a (e k ) = -V^ c 



£fc ~ AV 



(15) 



(16) 



Here /j£ c is the spin-dependent electrochemical potential 
defined by 



Ha 



eE ■ r + [i a 



(17) 



and the chemical potential fi a consists of terms in the 
zeroth and first orders of the driving forces: 



t'o 



(18) 



The collision term is written aJ^ 



dfa 

dt 



fe' 



[-W kk .„f a (k) + WW/ CT (fc')] , (19) 



where W^k'a is the rate of transition from ka to k'a and 
has the contribution from the normal scattering, W kk , a , 
and that from the skew scattering, W kk , a : 



Wi 



fcfc'r 



w. 



kk r < 



(20) 
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with 



We seek the solution for of the form 



w. 



kk'< 



W n (£ k ,e)5(e k ,-e k + eE-Ar), 

a sine W ss {e k , 9)6(e k , - e k + eE ■ Ar), 



(21) 



Here 9 is the angle of fc' relative to that of fc. Since we 
retain only terms up to 0(a), W ss (s k , 9) representing the 
skew scattering is O(a), while W n (e k , 9) due to the nor- 
mal scattering has no dependence on a. Both W n (e k ,9) 
and W ss (e k , 9) are an even function of 9. The delta func- 
tion expresses the conservation of energy, in which we 
take into account the potential energy shift due to the 
position change in the side jump at the scattering from 
ka to fc'c, 



Ar = —2aah(k' — fc) x e z 



(22) 



where e z = (0, 0, 1) and the vector fc should be regarded 
as a three-dimensional vector with vanishing z compo- 
nent, fc = (k x ,k y ,0). Note that the functions W n (e k ,9) 
and W ss (e k , 9) are defined in the absence of E where the 
difference between e k and e k ' is absent. 

The collision term is separated into four components, 



dt 



— CnO + Cnl + CssO + C 



ssl ; 



(23) 



with 



(24) 



C n o=£^'a [/ (0) (^)-/ (0) (e fe ) 

k' 
k' 

C ss0 = J2w s k s k , a [-/ (0) (e fe <)-/ (0) ( £fc ) 

k' 



We retain terms up to 0(E) and those up to 0(a). Then 
we immediately have C ss o = since the side jump Ar giv- 
ing terms of 0(a 2 ) in C ss o is to be neglected and the inte- 
grand of Csso becomes an odd function of 9. On the other 
hand, C n o is not zero in the presence of the side jump. 
The side jump gives the difference between f^(e k >) and 
f^{s k ) of C n o in two ways. One is from the difference 
in the kinetic energy s k , which comes from the potential 
energy shift and the energy conservation at the scatter- 
ing. The other is from the difference in the distribution 
between two points separated by Ar, which is described 
in the local equilibrium by the difference in and that 
in T e . Such considerations give 

/ (0) M - f {0) (s k ) = ^(e kl s k ) + ^Vav • Ar 



de k 

a/(°) 

dT c 



VT C • Ar, 



(25) 



and s k i — e k = —eE ■ Ar using the energy conservation. 



fP(k) = - d f^hk-V a (e k ), 
oe k 

and substitute this form into C n i and C ss i. 
straightforward calculation gives, for e = e k , 



(26) 
Then a 



dt 



de 



hk 



aV a {e) x e z 2aae z x F a {e) 



.(27) 



T n {e) Tss(e) r n (e) 

The first and second terms in the square brackets come 
from C n i (the normal scattering) and C ss i (the skew scat- 
tering), respectively, with r n and r ss defined by 

— -J2 5 ( £ k' -e)Wn(M)(l-cos0 



1 

T ss (s) 



k' 



8{e k , -e)W ss (e,9)sin 2 9. 



(28) 
(29) 



Note that r ss (e) can be negative since Wss(e,#) starts 
from the third order in the expansion with respect to the 
impurity potential.^The third term comes from C n o (the 
side jump) and is induced by the gradient of the chemical 
potential and that of the electron temperature as well as 
the electric field. Substituting the drift term Eq.( 15 ) and 



the collision term Eq.(27l into the Boltzmann equation 



Eq.(ll) gives the following equation for V cr (e): 
F CT (e) _ V a {e) aV a (e) x e z 2aae z x F a (e) 



T SS (£) 



Tn(£) 



(30) 



Up to the first order of the spin-orbit coupling constant, 
a, V a (e) is obtained to be 



V a {e)-- 



Tn{e) 



m 

-2aaF a (e) x e z 



F a (e)~a^\F a (e)xe z 



(31) 

Substituting this formula of V„{e) into that of in 
Eq.(26), we obtain the distribution function, / CT (r, fc), in 
Eq.(12) in the presence of the electric field, the chemical 



potential gradient, and the temperature gradient. 

C. Current densities and transport coefficients 

The number current density of spin-cr electrons is de- 
fined by 



(32) 



where the summation is taken over spin-cr electrons in the 
area S, and «j is the velocity operator of the ith electron 
given by 



Vi = — + 2aaVv 2 v{ri) x e z 
m 



(33) 
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The second term of V{ comes from the spin-orbit interac- 
tion induced by the potential due to the electric field and 
impurities, V2u( r i), and reduces to —2ao{dpi/dt) x e z in 
0(a). The brackets in Eq.(32) take the average with 
respect to the wave packet in the steady state. In the 
steady state the acceleration by the electric field is bal- 
anced with the deceleration by the impurity potential 
when each wave packet travels through the system, that 
is (dpi/dt) av — 0, which leads to the vanishing contribu- 
tion from the second term of Vi to the current. This semi- 
classical argument made by Hankiewicz and Vignale^ 
has been supported in terms of a rigorous density-matrix 
formalism by Culcer et alW^ The first term of Vi gives 



3 = 



hk 



s ^ 

k 



e k) 



(34) 



Here the contribution from vanishes since de- 
pends only on the magn itude of k. Substituting the ex- 
pression of /< 



(i) 



Eq.( 26 ), we have 
r = ( P eV a (e)) c 



(35) 



where p is the constant density of states per unit area 
per spin for two-dimensional electrons and the brackets 
represent the statistical average for spin-cr electrons: 



(••■>. 



de ■ 



o 



d/W( £ ,/v,r c ) 

de 



(36) 



The heat current density is obtained in a similar manner 
as 



(peV<r(e)(e - /v))<r- 



(37) 



In the linear-response regime, each component of the 
number current density j na is a linear function of com- 
ponents of thermodynamic forces, and the same is the 
case for j qa . The thermodynamic force corresponding to 
each current density is obta ined from the expression of 
the entropy productionj^^to be — T" 1 V/j£ c for j nrT and 
-T e " 2 Vr c for j qrT . Therefore the linear relations between 
the current densities and the thermodynamic forces are 
written as 



with the transport coefficients 



-V/4 C \ 



XX 



xy 

yy 



* = 1,2, j = 1,2. 



(38) 



(39) 



The common factor T^ 1 of the thermodynamic forces is 
absorbed in the transport coefficients. Since the 2DES 
in our model is isotropic in the xy plane, the transport 
coefficients have the following symmetry relation: L l J° — 



and L«" = 



yx 



The expression for each transport coefficient is ob- 
tained by substituting the formula of V a (e) in terms of 



the thermodynamic forces, Eq.(31) with Eq.(16), into 
those of the current densities, Eqs.(35) and (37 1. The 
obtained expression is 



L ija 



(40) 



with = x,y and v = x,y. Here L fiU (e) is the contribu- 
tion to the conductivity from electrons having energy e. 
Diagonal components 



L xx {e) = Lyy(e) = pe 



T-n(e) 



(41) 



have the form of the Drude conductivity divided by e 2 , 
while off-diagonal components 



L xy (e) = -L yx (s) = pea 



niT ss (e) 



2a 



(42) 



are due to the spin-orbit interaction. The first term in 
the square brackets is the contribution from the skew 
scattering, while the second term is that from the side 
jump. Note that the spin-orbit interaction gives rise to 
all off-diagonal transport coefficients in L lla , L 12a , L 21a , 
and L 22lJ . 

When the 2DES is degenerate (/v » k^T e ), 



L 



L^uiPa); L 1 



-(k B T c f 



dL l _ lu (e) 
de 



e=fi„ 



(43) 

Therefore the Mott relational holds, that is, the ther- 
moelectric conductivity tensor, L 12 ° , is proportional to 
the energy derivative of the electric conductivity ten- 
sor, Lj^J 7 . In addition, the diagonal electric conduc- 
tivity reduces to the Drude conductivity, e 2 L^ CT = 
n a e 2 T n (iJ l<T )/m, where n a is the density of spin-cr elec- 
trons. 

In the discussion of the spin Nernst effect as well as 
the spin Hall effect, it is convenient to reorganize the 
number and heat current densities for both spins into 
the spin current density, j s , the number current density 
ji n , and the heat current density, j q , as follows, 



j a =j nt +3 ni , 
m — -,-qt i »-q4- 



(44) 



where we have used the notation a =^, \, instead o f a = 
+1,-1. The corresponding thermodynamic forceiP^^ 
are — T e -1 V/4 c , -T^V/^, and -T- 2 VT C , respectively, 
with 



a4 = (a4 c + /4 c )A 

The linear relations now become 

L ss L sn L s q \ / _V/4 C 

L ns L nn £nq -V^ c 

L qs L qn L qq J \-T- l VT 









-( 







(45) 



(46) 
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where 

L iit +L m L m_ L m L 2it_ L au 

L ss = , L sn = , L sq = , 

4 2 2 

L ns = L sn , L nn =L n ^+L in , L nq = L 21t +L 2U , (47) 

L qs = L sq , L qn = L nq , J L qq = L 22 ^+ J L 22 ^. 

In the following we employ the condition satisfied in 
nonmagnetic systems, that is, the chemical potentials for 
both spins are the same in equilibrium, = fiy' . Then 
we have L l J^ = Ujj and L % J£ 



L l J}. With use of these 



relations, we confirm that the Onsager relatiorP^Ml is 
satisfied, that is L^J — LPJz- In addition we have 



L iJt , L ijl = 2Lt j h, I 



L m _ £«■!■ _ 2L*J£j, J 



1 
I 

I 

— I 



(48) 



from which we find that L SD (~ L DS ) and L sq (= L qs ) are 
proportional to J, while the other matrices in Eq.(46) 
are proportional to I. Therefore we can separate Eq.(46) 
representing the linear relations into the following two 
equations: 




and 









-( 







rss rsn rsq N 
xx xy xy 
r ns r nn j nq 
yx ^yy n yy 
L qs L qn £qq 
yx ^yy ^yy/ 



j ss j sn T sq 
yy yx yx 

j ns r nn r nq 
xy xx xx 

r qs rqn r qq 
ie re '-'xx/ 



(49) 



(50) 



These equations indicate that, in nonmagnetic systems, 
the spin current (say, along the a; axis) is coupled only 
to the perpendicular component of the number and heat 
currents (along the y axis). 



III. SPIN NERNST EFFECT 

A. Calculation of the spin Nernst coefficient 

We consider a state in which all current densities are 
uniform in a rectangular sample in the xy plane. In this 
state the thermodynamic forces are also uniform as de- 
rived from Eqs.(49) and (|50|). We apply a uniform tem- 



perature gradient along the x axis (V X T C = const. =^ 0, 
V y T c — 0), under the condition that both the number 
current and the spin current are vanishing (ji s = and 
j n = 0). The spin Nernst effect in the absence of the spin 
relaxation is the appearance of a uniform gradient along 
the y axis of the chemical-potential difference between 
up and down spins, /x| c = fj^ — /zj., proportional to the 
applied temperature gradient along the x axis: 



N S V X T C 



(51) 



Here we call N s the spin Nernst coefficient. 

To obtain the formula of iV s in terms of transport co- 
efficients, we write the conditions of if, — and j° = 
in terms of the thermodynamic forces using Eq. ( 50 1 and 
eliminate V, 



„. Then we obtain 
1 L\ 



N s = -, 



" nn t sq 
J xx yx 

'J 1 T nn T ss 

e -^xx^yy 



L sn L nq 

yx xx 

rsn [ ns : 

yx xy 



which becomes, in the first order of a, 



9 rllfr21t 

j\ j ^ XX XV 

N s = 

T c It lit 



L llf£21t 

xy xx 

2 ■ 



(52) 



(53) 



On the other hand, Eq.(49) with j% = 0, j° = 0, and 
\7 y T c = gives V^cc = V Vj,/i* = 0, and $ = 0. In 
particular \7 x /i s cc — means that no spin accumulation 
is generated in the same direction as the applied temper- 
ature gradient. 

In calculating the spin Nernst coefficient, we consider 
the degenerate electron gas in which the equilibrium 



chemical potential, fi = /ij 



(0) 



(o) 



is much larger than 



fceTo. In this case, using Eq.(43), we obtain 



27T 2 HT e 

A s = -2-2 



r ss (/i) r BB {n) 



with 



and 



r n (fj,) 2ma 



dT n (e) 



ds 



(54) 



(55) 



(56) 



e=fi 



The first term of n s in Eq.(55 1 comes from the skew scat 



tering, while the second term is from the side jump. Han- 
kiewicz and Vignale^l have shown, in the calculation for 
a model impurity potential, that the energy dependence 
of r ss is smaller than that of r n at the Fermi energy. 
Therefore the term with n s is dominant in N s in Eq. ( |54[ ) . 
They have also shown that r ss is negative (positive) for 
repulsive (attractive) impurity potentials.^ On the other 
hand, the spin-orbit coupling constant, a, is positive for 
semiconductors. 

According to Kohn and LuttingerpB both l/r n and 
l/r ss are proportional to the impurity density, n mip , up 
to the third order in the expansion with respect to the 
strength of the impurity potential. Therefore we employ 
this proportionality by considering weak impurity poten- 
tials. Then in Eq.(55 1 the first term of n s from the skew 
scattering is independent of 7ii mp , while the second term 
from the side jump is linear in ni mp and the coefficient is 
negative. In the case of the repulsive impurity potential 
where the first term of n s is positive, n s changes its sign 
when ni m p is increased, as shown in FigjTJ while, for the 
attractive impurity potential, n s is negative at any value 
of n imn . 



7 



(a) repulsive impurity potential 







h 



l imp 



tl. 



(b) attractive impurity potential 







h 



l imp 



FIG. 1. Normalized spin Nernst coefficient, n s in Eq. ( 55 1 and 
the normalized spin Hall coefficient, h s in Eq.(59l 
tion of the impurity density, ni mp , for (a) 



as a func- 



a repulsive impurity 
potential and for (b) an attractive impurity potential. 



B. Comparison with the spin Hall coefficient 

We compare the dependence of the spin Nernst coeffi- 
cient on the impurity density with that of the spin Hall 
coefficient derived in ReffTBI We apply the number cur- 
rent along the x axis, while we keep the spin current 
vanishing. We also set the condition that the current 
along the y axis is vanishing for both the number and 
the spin, and that the electron temperature is uniform. 
The condition of j s y — with V X T C = in Eq.(50l gives 
immediately 



V,Xc = ff s V x /C, (57) 
with the spin Hall coefficient H s given by 



' L ss - z r llf 
yy ^xx 



When the electron gas is degenerate, we obtain 

r n (/i) 2ma 



H s — 2h s , /i s 



(58) 



(59) 



which reproduces the dependence of the spin Hall con- 
ductivity on r n and r ss derived in 2DES by Hankiewicz 



and Vignale.^ Comparing this formula of h s with that of 
n s in Eq.(55l, the difference appears only at the sign of 



the second term from the side jump: the slope of h s as a 



function of 



''imp 



is positive, while that of n s is negative. 



Therefore the change in sign of h s with ni mp appears for 
the attractive impurity potential as shown in FigjT] 

The sample in the experiment by Sih et al.Pin which 
Si donors are doped in the quantum well, corresponds 
to the attractive impurity potential in FigfTJb). Accord- 
ing to the calculation for this sample by Hankiewicz and 
Vignalep^ the contribution to the spin Hall conductivity 
from the side jump is comparable in size to that from 
the skew scattering. Therefore we expect that the sign 
change of the spin Hall coefficient should be observed if 
the density of Si impurities in the well is changed around 
the value of the Si density used in the experiment. 

A Be impurity in GaAs is known to act as an accep- 
tor. Therefore doping Be in a quantum well introduces 
the repulsive impurity potential for the 2DESp^ In this 
case it is expected, according to the calculated result in 
FigjTJa) , that the spin Nernst coefficient changes its sign 
as a function of the density of Be impurities. 



IV. CONCLUSIONS 



We have studied theoretically the spin Nernst effect 
due to the spin-orbit interaction in the extrinsic origin in 
two-dimensional electron systems (2DES) in the xy plane 
by employing the Boltzmann equation. We consider a 
2DES confined within a symmetric quantum well with 
delta doping at the center of the well. We have shown 
in such a 2DES that the spin-orbit interaction, including 
that induced by the impurity potential, is diagonal in 
the z component of spin because of the symmetry of the 
system. 

In this model of 2DES we have investigated the de- 
pendence of the spin Nernst coefficient on the sign of 
the impurity potential and on the impurity density, and 
compared the result with that of the spin Hall coeffi- 
cient. We have found that the spin Nernst coefficient 
changes its sign as a function of the impurity density in 
the case of the repulsive impurity potential, while no sign 
change occurs for the attractive impurity potential. On 
the other hand, the spin Hall coefficient changes its sign 
in the case of the attractive impurity potential, as shown 
already by Hankiewicz and VignaleP^ The sign change 
of each coefficient occurs due to the cancellation between 
the skew-scattering contribution and the side-jump con- 
tribution. 
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